Mesoscopic conductance fluctuations in graphene samples 
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Mesoscopic conductance fluctuations in graphene samples at energies not very close to the Dirac 
point are studied analytically. We demonstrate that the conductance variance ([SG] 2 ) is very sensi- 
tive to the elastic scattering breaking the valley symmetry. In the absence of such scattering (disorder 
potential smooth at atomic scales, trigonal warping negligible), the variance ([8G] 2 } = 4{[5G] 2 ) mcta i 
is four times greater than that in conventional metals, which is due to the two- fold valley degeneracy. 
In the absence of intervalley scattering, but for strong intravalley scattering and/or strong warping 
([SG] 2 } — 2([(5G] 2 ) mc tai- Only in the limit of strong intervalley scattering {[<5G] 2 ) = ([<5G] 2 ) mc tai. Our 
theory explains recent numerical results and can be used for comparison with existing experiments. 

PACS numbers: 73.63.-b, 72.15.Rn, 81.05.Uw 
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Introduction. Graphene (a monolayer of graphite) is a 
novel material P, 0, 1, 0| with the Dirac electronic spec- 
trum. A lot of progress in the theoretical understanding 
of clean graphene has been made so far (see e.g. Ref. [H). 
For many interesting effects, however, disorder plays a 
significant role. A peculiar feature of disordered graphene 
is that its physical properties are sensitive to the scatter- 
ing processes breaking the valley symmetry. 

This sensitivity has been revealed in the behavior of 
the weak localization (WL) correction to conductivity 
[SUBS EH- Another famous phenomenon due to disor- 
der that goes along with WL are the mesoscopic conduc- 
tance fluctuations (CF). CF with variance ~ e 2 //i are ob- 



served in graphene samples experimentally 0, 11 , 12, R 
although a detailed analysis has not been reported. 

Numerical investigation of CF in graphene has been 
undertaken recently in Ref. [3] ■ The authors have rather 
unexpectedly found that CF in graphene were consider- 
ably stronger than those in conventional metals IH, 16| 



and the variance did not seem to be universal. No clear 



explanation of this effect was given in Ref. |14j , although 
it was argued that the unusual behavior might be due to 
percolation effects. 

Here we develop an analytical theory of conductance 
fluctuations in diffusive graphene samples at energies not 
very close to the Dirac point. The results we obtain ex- 
plain the findings of Ref. [3] and can be directly used 
for comparison with the experiments. 

Model. We consider a general microscopic model of 
disorder in graphene. The single-particle Hamiltonian of 
the system is (we put h — 1 and recover it later on) 



H r = H + H w + V(r), H = -iv(l 



KK' 



_AB 



)dx- (l) 



Here H w = fi w T^ K ' ® [r x (d^-d 2 )-2T y d x d y ] AB describes 
weak trigonal warping and A = x, y. The Hamiltonian H r 
is a matrix in the tensor product KK' ® AB of the valley 
(KK') and sub-lattice (AB) spaces and 1, t XiViZ are the 
unity and Pauli matrices. The random disorder potential 



V(r) is Gaussian with the correlation function 
(V-(r) <g> V(r')) = {A 1 ® 1 + A?T W ® T kl }S(r - r'), (2) 



where 1 = l KK ' <g> 1 AB , T kt = r 



KK' 



, k, I — x, y, z. 



For a given Fermi energy e the scattering rates (inverse 
scattering times) are defined as 

(7o,7**,7*-L,7-U,7-L_l) = nv(A Q ,Al, K z x yl Af y , A^^), 

where v = e/(2irv 2 ) is the density of states per one valley 
and one spin. 

In Eq. ([2|), the term <x 70 arises from remote charge 
impurities in the substrate, the field of which varies 
smoothly at atomic scales, while the rest of the terms 
describe various atomically-sharp defects that break 
the valley symmetry. Being diagonal in ATAT'-space 
(oc rf K ), the terms oc j zz and oc -f z ± do not involve the 
intervalley scattering, but do lift the valley degeneracy 
by acting differently on the valleys. Such terms describe 
the intravalley scattering, whereas the terms tx ^y± z and 
ex 7^ are due to 

Calculations. We calculate the correlation function 

F a0 , yS (Ae,H,AH) - (6a a p(e + Ae, H + AH)6a y5 (e, H)) 

(3) 

of the conductivities a a p(e + Ae, H + AH) and cr 7 5(e, H) 
taken at the Fermi energies e -I- Ae, e [rjj and magnetic 
fields H + AH, H (a, f3,j,S — x,y and 8a = a — (a) ) . 

We use the averaging technique developed for conven- 
tional disordered metals [l8[ . We assume (i) weak disor- 
der, A /v 2 <C 1, and (ii) diffusive regime, i.e., that the 
mean free path I — u/70 is much smaller than the size of 
the sample and the valley-symmetric rate 7o_js dominant, 
70 ^> 7{z,±}.{2,±}- As it was shown in Ref. [9( one should 
first renormalize the velocity v and constants A [Eq. |(5J)] 
solving renormalization group equations and then use 
them for calculating the localization corrections. This 
procedure is valid so long as e > eo exp(— ttv 2 /Aq), where 
eo is an atomic-scale energy. The same can be done for 
the correlation function (|3|) and we further assume that 
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FIG. 1: Diagrams for the conductivity correlation function 
F a p n g(Ae, H, AH) [Eq. Gray stripes denote diffusons 

and Cooperons, rendered with lines blocks are Hikami boxes, 
see Fig. The diagrams (c),(d) with the substitution a *-* /3, 
7 <-» 5, i? . <-> A must also be considered. 
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FIG. 2: (a),(b) Diagrammatic representation of the inte- 
gral equations for the diffuson and Cooperon. (c),(d) Hikami 
boxes. The current vertex renormalized by disorder (dark 
triangle) equals j a = 2ev(l KK ® t^ B )- 



v and A's have been renormalized. Under these assump- 
tions, the calculations for graphene generalize those for 
ordinary metals 15, 16| and F a t7 s(Ae, H, AH) is given 



by the diagrams in Fig. [TJ 

The arising disorder-averaged products of the exact 
retarded (R) and advanced (A) Green's functions Q R ' A 
yield the diffusons and Cooperons defined as: 

t) u {T,i>) = (g R (e + uj,v,v')®g A (£y,v)), (4) 

C u (r, r') ee (g R (e + u, r, r') ® G A (e, r, r')). (5) 

They satisfy the integral equations represented diagram- 
matically in Fig.[2[a),(b). These equations possess a non- 
trivial matrix structure acquired from the correlator ([2|) 
(dashed lines) and disorder-averaged (Q R ' A ) (fermionic 
lines). Consequently, there exist "high-energy" modes of 
D and C with gaps ~ 70 as well as "low-energy" modes, 
which gaps do not contain 70 @, @] • Giving much greater 
contribution in the diffusive regime, only the latter low- 
energy modes are of interest here. 

The effect of trigonal warping can be taken into ac- 
count in Q R ' A up to the second order in H w . As a result, 
the self-energies of the Green's functions {Q R ' A ) acquire 
the warping rate 7 W = (fi w e/v 2 ) 2 /7o and the dashed line 
in Fig. [2a), (b), in addition to the correlator ([2]), also 
represents the "warping term" 



7- 



D,C 



±7w(t. 



KK' 



1 



AB\ , 



ft 



KK' 



1 



AB\ 



with + and — for the diffuson [£),(a)] and Cooperon 
[C,(b)], respectively. It appears that, in the low-energy 
diffuson/ Cooperon subspaces, the matrix structure of the 
warping term 7^' C is identical to that of the intravalley 
scattering of type j zz = j zz T zz ® T zz in Eq. @. There- 
fore, the effect of warping on the low-energy diffusion 
modes of D and C is not any different from that of the in- 
travalley scattering cx 7 22 and the effect of warping could 
be taken into account by the substitution 7 22 — > 7 22 +7 w . 

Resolving the matrix-structure of the equations in 
Fig- Ufa), (b), we obtain: 



D u (t, r') = — { [D° + Dl + 2D 2 f] 1 <g> 1 + [£>° + - 2D 2 f] r z (g>r z + [£>° - £> A ® r A } <g> (1 <g> 1 + r k ® r k ) (6) 
a(r, r') = ^ { [2C 2 / 3 + Cl + C£] 1 ® 1 + [2C 2 / 3 - C£ - C£] r 2 <g) r 2 + [C* - C°]r A ® r A } (1 ® 1 - r fc ® r fe ) (7) 

I 



In Eqs. © and (0, the tensor products are ordered as 
(i? (g) Ajifjc ® {R® A)ab, A = y, and k = x,y, z. 
The diffuson/Cooperon components = £>^(r, r') and 
C£, = C^(r, r'), i = 0, 1,2,3, satisfy the equations 

{-z^-^V^ c +r. 1+ 7 lncl }(^,C)L(r,r') = <5(r-r'), (8) 

where V D = V-i(e/c)AA(r), Vc = V-i(e/c)[2A(r) + 
AA(r)], the vector potentials A(r), AA(r) correspond 
to H and AH, respectively, D = v 2 /-fo is the diffusion 



coefficient, and 7i no i is the inelastic scattering rate due 
to, e.g., electron-electron or electron-phonon interactions. 
The elastic scattering rates Ti due to disorder equal 

r = 0, ri = 4 7± , r 2 = r 3 = 2 7± + 2 7z + 2 7w , (9) 

where the total intervalley 7j_ = 7± 2 + 27j_j_ and intraval- 
ley 7z = 7z 2 + 2j z ± scattering rates were introduced. 

The Cooperon in the form of Eq. ([7]) has been obtained 
earlier 0, @| , whereas the form © of the diffuson is ob- 
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tained here for the first time. Note that for a given i the 
rates [Eq. |9])] entering the corresponding diffuson 
and Cooperon modes are identical. The fact that the 
insensitive to various phase-breaking phenomena diffuson 
D does contain the rates j z , 7j_ , and 7 W (lij means that 



the effects of the intravalley and intervalley scattering 
and of the trigonal warping should never be understood 
as a suppression of electron interference alone. 

Results. Calculating the diagrams in Fig. [TJ for the 
correlation function of conductivities [Eq. ([3])] we obtain 



F a ^ s (Ae, H, AH) = (Sa aP (e + Ae,H + AH)Sa j5 (e, H)) = (2. 2 v e 2 D) 2 J ^ 



deals' dn(e) dn(e') 1 



de de> S 2 



J drdr ' E {s cn 5p S \Dt(Ty)\ 2 +5 aB 5fr\Ci(T,r')\ 2 + ~5 a/3 6 lS Re[Di(r, r')D'l{r' , r) + C£(r, r')^(r', r)]| , (10) 



r 



where u> = e — e' + Ac, S is the sample area, n(e) — 
l/[exp(e/T) + 1] is the Fermi distribution function, and 
the factors 2 S and 2 V originate from the dimensionality of 
the spin and valley spaces (the indices s and v emphasize 
their origin). 

Equation (fl~0|) , together with Eqs. j6])-([9]), constitutes 
the main result of our work. The key feature character- 
izing graphene is that different diffuson and Cooperon 
modes i enter Eq. (JTUJ) . The magnitude of mesoscopic 
fluctuations is thus determined by the strength of elastic 
scattering processes breaking the valley symmetry. When 
all such effects are negligible, the result (|T0|) for graphene 
is (2„) 2 = 4 times greater than that for conventional met- 
als [H, [3| due to the two-fold valley degeneracy. Note 
that at H = AH = one has D l u = C l u and for a given 
i the diffuson and Cooperon contribute equally. 

To be specific, below we consider the case of a rect- 
angular sample with length C x and width C y , occupying 
the area < x < C x , < y < C y , and attached to 
ideal leads at x = and x = C x . The conductance 
G = G xx in the x direction is related to the conductivity 
a xx as G = a xx C y /C x . From Eq. ((TUJl, at T, 7 inc i < e*, 
where e* = ir 2 D/C 2 is the Thouless energy for the 
x dimension, for the conductance correlation function 
T{Ae) = (SG(e + Ae)5G(e)) = {C y /C x ) 2 F xx>xx (Ae,H,0) 
we obtain 



F(Ae) = a H 



2^ 

2^ 



1 



x EE { 2 I^X e (q)| 2 + Re{D* Ae ( q )} 2 } ,(H) 

i=0 q 

where D l Ae (q) = l/(— iAe + TJq 2 + Ti) are the spatial 
eigenmodes of the diffuson, q 2 = q 2 



= 1,2, 



and 



<ly 



<,yj Cy, 



n y = 0, 1, ... . In 



Eq. (fTTj) . the factor an accounts for the sensitivity of 
the Cooperons to the magnetic field in the two limit- 
ing cases: an — 1 for H <C H* and an — 1/2 for 
H ^> H* , where H* = (c/e)/C 2 . The conductance vari- 



ance F(Ae = 0) = ([<5G] 2 ) following from Eq. (JTT]) equals 

n 2 3 



\SG} 2 ) = 3a 



H 



2 S 2 v e 



2nh 



J2K(jCi,jC x ,jC y ){12) 



oo oo 

"74 ^2 



7T 4 £t 



n x — 1 n v — 



l 

C 2 



C X Cy 



GL3) 



where C 2 = Ti 2 D/Ti. For both narrow [C y <C C x ) and 
wide {C y 3> C x ) samples, the contribution of a given 
mode i to Eq. (fT5|) is unsuppresscd if Ci ^> and equals: 



72.0 = 7£(oo, C X ,C V ) = 



1/90, 

C(3)V(47r 3 £ x ), £ 



ZCy aCj: 



Wide samples are thus more attractive for the observa- 
tion of unsuppressed CF. In this case the length Ci has 
to be greater than only the shorter dimension C x (equiv- 
alently, Ti <C e*), but can be arbitrary compared to C y . 

The limiting cases of Eq. (fT2|) for different strengths of 
the scattering processes can be summarized as follows: 



([SGf 



'-aJ\8G\ 



{[5Gf 



12a H 



2nh 



where the conductance variance ([5G] 2 ) m for a conven- 
tional metal TH, 16| was introduced. The coefficient a 7 
gives the number of diffusion modes i that contribute 
(i.e., for which r i < e*) to CF, see Table HI As follows 
from Eq. ©, the mode i — ("pseudo-spin singlet") 
is unaffected by any of the scattering mechanisms, the 
mode i = 1 ("triplet, 0") can be suppressed by the inter- 
valley scattering only, and the modes i = 2, 3 ("triplet, 
±1") can be suppressed by both intervalley and intraval- 
ley scattering and by trigonal warping. Note that trigo- 
nal warping does affect CF, in the same way as intravalley 
scattering does. 

(i) When all the effects are negligible, r y± 7 -f z ,lw <C e*, 
all four modes contribute equally, a 7 = 4, and ([5G] 2 ) = 



4 





7z -C and ~/„ *C ej 


7z > ej or 7™ > ej 




4 


2 


7-L > e£ 


1 


1 



ing to the conductance variance in graphene for different in- 
tervalley 7x, intravalley 7 Z and trigonal warping 7 W scatter- 
ing rates. The value of q 7 also gives the ratio of the con- 
ductance variance in graphene to that in conventional metal, 

([(5G] 2 ) grapheme = Q 7 ( [5G] 2 }metal , See Eq. JUJ . 



4([(5G] 2 ) m is four times greater than that for a conven- 
tional metal. This is explained by an additional two-fold 
valley degeneracy described by the factor 2 V in Eq. (fl"4l) . 

(ii) If the intervalley scattering is weak, 71 < e*, but 
either the intravalley scattering or the trigonal warping 
are sufficiently strong, j z ^> e* or 7 W ^> e*, then the two 
modes i = 0, 1 contribute, while the modes i = 2, 3 are 
suppressed. In this case a 7 = 2 and ([<5G] 2 ) = 2([<5G] 2 ) m 
is two times greater than that for a metal. 

(iii) Finally, if the intervalley scattering is strong, 
7_i_ <C e*, and the intravalley scattering ^ z and trigo- 
nal warping 7 W rates are arbitrary compared to e*, then 
all triplet modes i = 1, 2, 3 are suppressed, and only the 
gapless mode i — contributes. In this case a 7 = 1 and 
([5G] 2 ) = ([<5G] 2 ) m coincides with that for a metal. 

The conductances fluctuations in graphene were stud- 
ied numerically in Ref. [14| (see Fig. 3 therein). For 
atomically-sharp disorder, a plateau ([^G] 2 ) p i atC au ~ 
([SG] 2 ) m in the dependence of ([<5G] 2 ) on the disorder 
strength Aq/v 2 was obtained, which clearly corresponds 
to the case (iii). For atomically-smooth disorder, a wide 
peak in the dependence of ([SG] 2 ) on Aq/v 2 with maxi- 
mum ([5G] 2 ) max ~ (4.5 - 5)([5G] 2 ) m at Aq/v 2 ~ 1 was 
obtained. We believe this situation corresponds to the 
case (i), the maximum value being close to our predic- 
tion. We emphasize that our theory, just like that of 
Refs. HI 3], requires both weak disorder {Aq/v 2 <C I) 
and the diffusive regime [I = w/70 <C C = min(£ a; , C y )}. 
The reason for having a pea k, rather than a plateau, for 
smooth disorder in Ref. 14] is that the range of Aq/v 2 , 
where both these conditions are met, is quite narrow. 
The diffusive regime is not reached until disorder becomes 
strong (Aq/v 2 > 1), while for smaller values of Aq/v 2 <C I 
the system is simply in the ballistic regime I > C. This 
is supported by direct check of parameters {1/ C ~ 0.1 for 
A /v 2 ~ 1 and thus l/C ~ 1 for Aq/v 2 ~ 0.1 ) and by 
an improving tendency (earlier upsurge of ([<5G] 2 ) with 
increasing Aq/v 2 <C 1) for larger samples (filled vs. open 
symbols). 

Our theory thus helps understand the findings of 
Ref. [IH without assuming the existence of percola- 
tion paths and nonergodicity as was done by the au- 
thors. The ergodicity implies equivalence of averag- 
ing over disorder and the Fermi energy (or magnetic 



field), (f[G(e)}) = UmAe^oo^/^e/a^'/IGCe')]. For 
f[G(e)] = G(e), this is clearly true, if F(Ae) as 
Ae — > 00, see Eq. ([Pf]) . This asymptotic of J-{Ae) is de- 
termined by the behavior of the diffusion modes D l Ae (q) 
at large energies Ae and, in this respect, graphene is not 
any different from an ordinary metal [2fJ. One can es- 



timate ^(Ae) oc J 



oc 



D 
Ac 



for Ae > e*,e*. 



The proof for f[G(e)} = [G(e)] m , m > 1, is analogous. 
Thus, the ergodic hypothesis for graphene holds. The 
violation of ergodicity in Ref. [3] occured near the An- 
derson metal-insulator transition, and might be due to 
the fact that averaging over energy was mixing extended 
and localized states. 

Conclusion. We have developed a theory of conduc- 
tance fluctuations in monolayer graphene samples. We 
expect our findings presented in the Results section to be 
also completely applicable to bilayer graphene samples. 

We thank SFB Transregio 12 for financial support. 
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